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INTRODUCTION 


The Tethered Satellite System (TSS) is designed to deploy a 
1.6 m diameter spherical satellite a distance of 20 km above the 
space shuttle orbiter on an insulated conducting tether. Because 
of the passage of the conducting tether through the earth's 
magnetic field, an emf is generated producing a positive 
satellite potential of about 5000 V. Electron flow under the 
influence of this high positive potential is the focus of the 
present analysis. 

The ionospheric parameters at TSS orbit altitude are; 
thermal velocity of electrons, 1.9 x 10 5 m/s, thermal velocity of 
the ions, 1.1 x 10 3 m/s, velocity of the satellite 8 x 10 3 m/s. 

The electrons, with a Debye length, A, D = 0.49 cm, spiral about 
the earth's magnetic field lines (0.4 Gauss) with a radius of 
about 3 cm and the ions spiral with a radius of 5 m. Under these 
conditions, the electron thermal energy, kT is 0.17 eV. The TSS 
satellite radius, r p is 163 Debye lengths. 

There is an extensive literature on the interaction of 
satellites with the near-earth ionospheric plasma. The space 
charge limitation to the electron current collected by a sphere 
at positive electrical potential was calculated by Langmuir and 
Blodgett (1924) . Parker and Murphy (1967) recognized the 
importance of the influence of the earth's magnetic field and 
used the guiding center approximation to calculate the electron 
current collected by a positive charged satellite. More recently 
Ma and Schunk (1989) have calculated the time dependent flow of 
electrons to a spherical satellite at positive potential 
utilizing numerical methods and Sheldon (1994) used similar 
methods to solve this problem for the steady state. 

In order to analyze some of the phenomena that occurred in 
the ionosphere during the TSS flights, it would be useful to have 
analytic expressions for these electron flows. The governing 
equations are very complex and an exact analytical solution is 
not likely. An approximate analytical solution is feasible 
however, and the results of one attempt are presented herein. 

MATHEMATICAL MODEL 

Electron flow to the spherical satellite is modelled here by 
the use of the cold plasma one-fluid momentum and continuity 
equations with the Poisson equation used previously [Sheldon 
(1994)]. In the present calculation, it is assumed that there 
will be a sheath region around the satellite devoid of ions due 
to positive satellite potential being much higher than the ram 
ion energy. Electrons drift toward this sheath along the earth's 
magnetic field lines with an average velocity (kT/m) 1/2 , 
neglecting the average initial velocity due to their ambient 
spirals. At the outer boundary of this sheath there will be a 
negative space charge potential barrier approximately equal in 
magnitude to the electron ambient thermal energy. The governing 
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equations are written in non-dimensional spherical coordinates 
(r,0 ,<j>) with the magnetic field aligned with the polar axis. The 
steady-state momentum equations for u r , u 0 and u., the r, 0, and 
<p components of the average velocity of the electron flow are 
then : ' 
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where r is in units of and the velocities are in units of u = 
(kT/m) , where m is the electron mass. The electric potential 
energy, V(r,0) is in units of kT and B = « /w = 0 . 2 , where u> is 
the electron cyclotron frequency in the earth^s magnetic field 
and o> e p is the electron plasma frequency. The continuity equation 
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where n(r,0) is the electron density in units of the ambient 
electron density, n A and Poisson's equation is 
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SOLUTION 

In order to get an approximate solution to Eqs(l)-( 5 ) the 


1 perturbation expansion in (l)-(5), 
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( 9 ) 


The equations to first order in B are 

[^( 1 , 6 ) u r ° ( r ) -F 1 (r, 0 ) ] =0 


u r °(r) J- [rufU'Q)] = 0 


( 10 ) 


Si^Mr,©) u* 1 ( r , 0 ) 

u°(r) [ — +— 5 -Bsin0] =0 

dr r 


(11) 


±J-[r 2 (nHr,Q) u z °(r) +n° (r) u z x (r,0) ) ] 

+ -^el [n ” (r>u ' 1(r ' e)sinei '° 


_d_ 

dr 


r r 2 _3vM_r L 0 )_ , 1 d 

dr sin0 60 


avMr,©) 

00 


sin0] =n 1 (r, 0) r 2 


(13) 


The solution to the zero order Eqs. (11) - (13) can be obtained 
from Langmuir and Blodgett (1924) . They considered a spherical 
electron emitter of radius r m outside of a collector with 
negligible initial velocity of the electrons leaving the emitter. 
The boundary conditions were V°(r m )=0 and (dV°/dr) rm =0. In the 
present case the outer boundary condition is taken at a virtual 
cathode of radius r 0 , where u r °(r 0 )=l and n°(r 0 )=l. 

With the boundary condition, u^ 1 (r o ,0)=O, Eq.(ll) gives, 

u 4> 1 (r,0)=B-^§(r o 2 -r 2 ) (14) 


The remaining first order perturbation equations are linear 
and allow the separation of variables. Defining u r 1= g r (r) g e (0) , 
u 0 1= h r (r) h e (0) , n 1= k r (r) k e (0) , V 1= f r (r) f e (0) , the separated 0- 
dependent Eqs. yield the following results 
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XXXIX-3 



(17) 
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P f (cos0) is the Legendre polynomial of the first kind. 

The r- dependent equations can be combined into a coupled 
set, of linear Eqs. which do not reduce to a standard form. An 
approximate solution is obtained by expanding the dependent 
variables in a Taylor Js series about r Q . Only the leading terms 
a ^" e use< ^ here. Combining with the 0— dependent functions, the 
final results largest order are 
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where y=ln(r m /r) and i= 47 tr 0 2 n A u A . 

The results for the azimuthal velocity, Eq. (14) , are 
compared with the numerical results of earlier work 
[Sheldon ( 1994 ) ] in Fig. 1. While the present method introduces 
considerable error in the actual values u., the radial profiles 
h^ve the correct shape. Similar results were obtained for 
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u r (r,0), V(r,0) and n(r,0) . 

CONCLUSIONS 

The perturbation method has allowed approximate 
determination of the electron flow in the proposed model. The 
previous numerical solution was more accurate, however the 
intention here was not accuracy, but a better understanding of 
the influence of the controlling parameters. This is available in 
Eqs(18) - (23) . 
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Figure 1. Radial profiles of the azimuthal electron velocity. 

V=1000. Nondimensional units 1 D and u A defined in text, 
(a) 0=22.5°, (b) 0=67.5°. 
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